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SUMMARY 

A method for use in the design of a complex engineering system by 
decomposing the problem into a set of smaller subproblems is presented. 
Coupling of the subproblems is preserved by means of the sensitivity 
derivatives of the subproblem solution to the inputs received from the 
system. The method allows to divide work among many people and computers. 

INHALTSANGABE 

Es wird vorgestellt, ein komplettes Engineering System zu entwickeln, das 
man in kleinere Untersysteme zerlegen kann. Das Zusammenlegen von 
Unter problemen kann erreicht werden mittels der Verf einerungen der 
Ableitungen von Unterproblemen , zu loesen empfangen durch Eingaben welche 
vom System werden. Die Methode teilt die Arbeit unter vielen Leuten und 
unter vielen Rechnern auf. 

INTRODUCTION 

Engineering systems are frequently composed of subsystems which are 
mutually coupled so that a modification in one affects the others and 
ultimately influences the performance of the whole. Designers of such 
systems are confronted with a closed loop situation whereby the decision 
making at the system level requires information about the subsystems that 
have not yet been designed and vice r ersa. The iterative approach 
commonly used to break the deadlock regies on past experience, judgment, 
and intuition as substitutes for the unavailable information and becomes 
increasingly inadequate as the systems grow in complexity and advance 
far away from the experience base. 
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The paper outlines a method for improving the iterative design of complex 
systems by making it more systematic and based on a set of coherent 
mathematical concepts. While rooted in mathematics, the approach 
deliberately avoids handing over the entire design process to the 
computer. On the contrary, the emphasis is on a broad work front of 
people and computers to combine the computer efficiency with the human 
intelligence indispensable in design. 

DECOMPOSITION 

The key to the propose 1 '* approach is a formalized decomposition of the 
large design problem into a set of smaller manageable subproblems coupled 
by means of the sensitivity data that measure the change of the 
subsystem design due to a change in the system design. Let ES be an 
engineering system composed of the subsystems SS^ , SS^, • • -SS^ . . .SS n as 

shown in Fig.1 (the abbreviations are defined in Table 1, and Table 2 
gives examples for the generic quantities in the context of aircraft 
design). The design variables are grouped in a vector SV for ES and the 
vectors DV. for SS . . The ES has a performance index PS that should be 

maximized within the system constraints collected in a vector GS. The ES 
imposes demands on each SS^. These demands are quantified by entries of 

a vector DS^ which depends on SV through analysis of ES. Suppose that 

each SS ^ is designed by manipulating DV^so that it meets its DS^, 

regarded as constants, while maximizing its safety margin SM.. 

representing a set of subsystem constraints GSS^ . These tasks separate 

for each SS^ can be carried out concurrently by whatever means the SS^ 

designers choose, including the appropriate analysis, optimization, and, 
also, judgment and experimentation. 

A new element required under the proposed approach is evaluation of the 
sensitivity of the maximum (optimum) SM^ to changes in in form of the 

optimum sensitivity derivatives SSM^/BDS^ At the ES level, these 

derivatives combined with the derivatives 3DS^/3SV in chain 

differentiation yield the sensitivity of SM^ to changes in SV in form of 

derivatives SSM^SSV. The maximum SM i and its derivatives show the ES 
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designer, with a linear extrapolation accuracy, how the change of SV that 
he controls will affect the SM i for each . Guided by this information 

and by the ES analysis, the ES designer can decide which variables in SV 
to change and by how much in order to move toward the goal of making all 

z 

the constraints GS and GS^ satisfied while maximizing the PS. The SV 

change will alter the DS^. Responding to that, the SS^ designers modify 

their designs and pass updated information to the ES designer who, then, 
changes the SV again, and so on. In this manner the ES and the SS i 

designers carry on a systematic iteration toward the improved system 
design, trading the data precisely defined in form of the DS^ , SNL , and 

their derivatives. Each designer works on a separate assignment with the 
control of PS vested in the ES designer while the SS^ designers focus on 

their SS i feasibility. The whole problem is decomposed yet remains 

coupled by the ES-SS i data exchange shown in Fig.1. 

OVERALL PROCEDURE 

Based on the above qualitative description, one may now formulate a 
step-by-step procedure to implement the decomposition approach. 

STEP 1. Initialize the system. 

STEP 2. Analyze the system. Calculate PS, GS, DS i , and 9DSV9SV. 

STEP 3- Design subsystems SS^ The DV., are manipulated within the upper 
and lower bounds, L^ and U^ f so as to find maximum SM. for given DS^ . The 
latter requires vector of equality constraints GE.. for those DS^ that are 
also functions of DV^ . These constraints enforce equality of the DS i 
values prescribed at the system level and computed as a function of DV i 
so that GE^-DS^ (SV)-DS^ (DV. )=0. Formally, the task may be formulated as 
an optimization 

max SM. (DV. ,DS. ) subject to constraints (1) 

DV. 1 

l 

GE.(DV.,DS.) = 0. L. S DV.S U 
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The output of the operation is: SMj.-fSI^) , and the optimal subsystem 
design variables, DV^ 

STEP 4. Analyze each SS^^ design for sensitivity to the inputs received 
from the system to obtain the 3SMV3DS^. 

STEP 5. Modify the SV to improve the system design. In this operation, 
one uses the 3DS./3SV, SM^ and SSM^/SDS^^ obtained in STEP 2, 3, and 4, 

to extrapolate each SM i as a function of the increment AV 

3SM . 3DS 

SMjUSV) • SM. . ^ 555- SSV (2) 

Improvement of the system design may be formalized as an optimization: 

a) max PS(SV) subject to constraints (3) 

SV 

b) GS(SV) £ 0, c) SM.(SV) £ 0 (for all i) 
d) L £ SV £ U 

in which the system level analysis provides the PS and GS, and the SM^ in 

eq.3c is approximated by eq.2. The bounds in eq.3d include "move limits" 
protecting the accuracy of the extrapolation in eq.3c. The above 
optimization problem may have no feasible solution within the move limits 
in eq.3d, if it begins with significant constraint violations in eq.3b 
and c. If a feasible solution can not be found, an acceptable outcome of 
eq.3 is a new design point moved as close to the constraint boundary as 
possible. The result of this step is a new SV defining a modified design 
of the system. 

STEP 6. Repeat from STEP 2 until all the constraints GS are satisfied, 
all safety margins SM i are non-negative, and the performance index PS has 

converged. 

In the above procedure, also shown in Fig. 2, the analyses in STEP 1 and 2 
are problem-dependent . The behavior sensitivity analysis required to 
obtain the 3DS i /3SV can be obtained by either a finite difference 

technique or, preferably, by a quasi-analytical method, e.g.,[1]. The 
optimizations defined by eq.1 and 3 can be carried out by any suitable 
algorithm capable to search an n-dimensional constrained design space, 
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e.g.,[2], although the use of a formal optimization method is a 
recommendation rather than a requirement. The optimum sensitivity- 
derivatives in STEP 4 can be calculated by means of the algorithms 
described in [ 3 ] , [M] , [ 5 ] , [6] . Extension of the above two-level algorithm 
to multilevel systems is given in [73 , and its application to aerospace 
systems is discussed in [8]. 

Initial tests including two-level [93 and three-level [103 structural 
optimizations showed satisfactory comparisons with the results obtained 
without decomposition. A status report on a multidisciplinary test case 
is documented in [11 3. It involves redesign of a wide-body transport 
aircraft wing for improved fuel consumption for a prescribed mission 
under constraints imposed by strength, aerodynamics, and aircraft 
performance requirements. Fig. 3 shows the aircraft and a three-level 
decomposition scheme devised for the problem. Research and development 
continue to learn more about the algorithm's convergence properties, 
sensitivity to the extrapolation errors and lack of synchronization among 
the subtasks, ability to handle direct couplings among the SS^s, ability 

to adjust to discrete or judgments.' decisions, and computer hardware 
dependence. 

SUMMARY AND CONCLUDING REMARKS 

The paper has presented a method for decomposing a large engineering 
design problem into a set of smaller subproblems. Each subproblem is 
self-contained so that all the subproblems can be worked on concurrently 
and a broad work front of people and computers can be developed. The 
method's testing to date has been satisfactory according to the 
references cited. Further development of the method is focused on the 
algorithmic details and continuing test applications in aerospace design. 
The development entails also the issues of the data management and 
utilization of parallel computation capabilit '-.-a being offered by modern 
computer technology. 
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Table 1 Summary of Generic Terms. 


DS. - vector of demand quantities imposed by the system on 

subsystem i . 

DV^ - vector of design variables for subsystem i. 

ES - (engineering) system. 

GE i - vector of equality constraints for subsystem i. 

GS ~ vector of system inequality constraints; an inequality 

constraint is defined as g=(DEMAND/CAPACITY)~1 , satisfied 
when g i 0, 

GSS^ - vector of inequality constraints for subsystem i. 

L, L. - vector of lower limits on SV, and DV. , respectively (move 

limits included). 

PS - performance index for ES (a scalar). 

SM^ - safety margin for SS^ (a scalar), defined as 

SM. =max(CAPACITY/DEMAND)-1 , or in terms of a set of g's 
(see GS): SM i «max(-g/ ( g+1 ) ) . 

SS^ - subsystem i. 

SV - vector of system design variables. 

U, U - vector of upper limits on SV , and DV , respectively (move 

limits included) 


Table 2 Examples of the Equivalents of the Generic Terms Typical 
for an Aircraft Application. 


DS - at the middle level: lift required of the wing; at the 

l 

bottom level: edge loads N , N , N on a wing cover panel. 

x y xy 

DV i - at the middle level: wing bending stiffness distribution; at 

the bottom level: detailed wing panel dimensions. 

E,S - aircraft, top (system) level; 

GE^ - at the middle level: wing structure weight prescribed at the 

top level; at the bottom level: panel spanwise membrane 
stiffness prescribed at the middle level. 

GS - runway length. 

GSS i - at the middle level: wing tip deflection; 

at the bottom level: panel local buckling. 

PS - fuel economy for a given mission. 

SS. - the wing box, middle level; the wing cover stiffened panels, 

third (bottom) level. 

SV - wing structural weight and airfoil thickness to chord ratio. 

3DS./3SV - derivative of wing lift with respect to structural weight. 

l 

3SPL/3DS. - derivative of wing panel safety margin with respect to edge 
loads . 
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